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We demonstrate how near-perfect entanglement (in fact arbitrarily close to maximal entanglement) can be 
generated between the end spins of an anti-ferromagnetic isotropic Heisenberg chain of length TV, starting from 
the ground state in the N/2 excitation subspace, by applying a magnetic field along a given direction, acting on 
a single spin only. Temporally optimal magnetic fields to generate a singlet pair between the two end spins of 
the chain are calculated for chains up to length 20 using optimal control theory. The optimal fields are shown to 
remain effective in various non-ideal situations including thermal fluctuations, magnetic field leakage, random 
system couplings and decoherence. Furthermore, the quality of the entanglement generated can be substantially 
improved by taking these imperfections into account in the optimization. In particular, the optimal pulse of a 
given thermal initial state is also optimal for any other initial thermal state with lower temperature. 
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I. INTRODUCTION 



Spin systems are important models to understand many- 
body physics and to investigate many interesting problems 
in the area of quantum information processing (TJ. Coherent 
control of such systems to achieve certain key tasks such as 
entanglement generation is one of the main problems in quan- 
tum mechanics. There have been several proposals already 
for entanglement generation between the end points of a spin 
chain both by distributing the entanglement through the chain 
112] O and entanglement generation through a quench followed 
by the natural non-equilibrium dynamics of the system |0}{6]. 
Unfortunately, the attainable entanglement in these schemes 
is not maximal and achieving an arbitrarily high entanglement 
between ending spins demands more control. Although there 
has been significant progress in both the theory and imple- 
mentation of coherent control, the microscopic nature of most 
quantum objects makes it difficult to control all elements of 
the system individually. Local addressing of individual spins, 
for example, remains a significant challenge. Applying fields 
that globally address the entire system is one way to avoid 
such problems |7|. On the other hand, if assuming locally ad- 
dressing ability of a many-body system, one interesting ques- 
tion is what is the minimal controlling resource that can per- 
form certain desired task. For example, universal quantum 
computation can be achieved through controlling a small part 
or gateway of the whole system [8 -12], when the system in- 
trinsic Hamiltonian provides the necessary non-local interac- 
tions for a successful control. Aside from not requiring access 
to the entire system, limiting the interaction of the controller 
with the system can be advantageous as quantum coherence 
is very fragile and any incoherent interaction with a macro- 
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scopic control apparatus can contribute to decoherence. Min- 
imal control is therefore likely to reduce controller-induced 
decoherence, which is a limiting factor for many quantum in- 
formation tasks. It has been shown that a wide class of cou- 
pled many-body systems can be manipulated as required by 
controlling a local subsystem |8|, and even for generic (not 
necessarily multi-partite) systems, control of a small subspace 
or even single transition is often sufficient to achieve full con- 
trollability [ 1 3 1 . In particular, the XXZ-Heisenberg model is 
controllable if one end spin can be fully controlled [9]. 

Controllability, however, only guarantees the existence of a 
control to complete a certain task, but does not tell us how to 
generate such a control, and even constructive controllability 
proofs do not usually result in efficient solutions. Efficiency 
is critical due to limited coherence time for many systems, 
which makes long pulse sequences problematic. Moreover, 
practical feasibility considerations impose various constraints 
on the permissible fields, from rise time to the complexity of 
its frequency spectrum. Consequently, control design has to 
be considered case by case and there is no general scheme 
that is optimal for an arbitrary system. Recently, there have 
been several proposals for achieving certain tasks in quantum 
information. For instance, an efficient transfer of quantum 
states can be achieved in an Ising chain, which is accessi- 
ble in Nuclear Magnetic Resonance (NMR), through global 
|7| and mixture of global-local [14] pulses. In |[T5l . a realis- 
tic scheme is proposed to achieve high fidelity state transfer 
by applying a sequence of two-qubit gates at one end of the 
chain and it was shown that the two qubit gates can be imple- 
mented by a proper sequence of switching the interaction of 
the two spins. Universal quantum computation along the XY 
spin chain can be achieved by controlling the two spins at one 
end of the chain |[T0lfT2ll . In JTD] one spin is controlled by two 
independent magnetic fields and a fast Hadamard gate operat- 
ing on the other spin. Alternatively, in |[T2l the two spins are 
controlled via three independent control fields, two of which 
are non-local interactions. Also, very recently, almost per- 
fect state transfer has been achieved by controlling the mag- 
netic field in the z direction over all spins in the ferromagnetic 
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Heisenberg chain |[T6l . 



In this work, we consider the problem of entanglement gen- 
eration between the two ends of an anti -ferromagnetic Heisen- 
berg spin chain when control is restricted to a uni-directional 
local magnetic field on one spin at either end of the chain. 
This type of control is not sufficient for full controllability of 
the iV-spin system lfT31 as both the system and control Hamil- 
tonian are excitation preserving, and thus the system decom- 
poses in into N + 1 excitation subspaces. For a Heisenberg 
chain, we shall demonstrate that with the Heisenberg inter- 
action and the local control Hamiltonian, the system is con- 
trollable in the particular subspace with the largest number 
of excitations, and thus if the target state is in the same sub- 
space as the initial state, it is reachable by such local con- 
trol. Through optimization, we find the time-varying con- 
trol field that generates an almost perfect singlet pair between 
two ends of the chain, after the system is initially prepared 
in the ground state. In general, calculations to generate the 
optimal control field for the largest excitation subspace of 
the anti-ferromagnetic chain are considerably more demand- 
ing than the more common single-excitation subspace calcu- 
lations 1 16 1, as the system state lies in the iV/2-excitation sub- 
space for N even, which has dimension Nl/[(N/2)l} 2 as op- 
posed to dimension N for the first excitation subspace, and 
has a complicated form being a superposition of many compu- 
tational basis states. Considering this scenario, however, has 
several benefits from a practical perspective: (i) most realiza- 
tions of spin chains are anti-ferromagnetic in their nature ifTTl ; 
(ii) in our entanglement generation strategy the initial state is 
the ground state of the system, which any system naturally 
takes when it cools down and thus there is no need for extra 
control to initialize the system; (iii) since the ground state of 
the system is highly entangled it can be considered as an ini- 
tial source for quantum information tasks. This helps simplify 
the control mechanism as multipartite entanglement naturally 
exists and the control simply has to convert it to utilizable bi- 
partite entanglement between the distant end spins. Despite 
the highly restrictive nature of the control and the dimension 
of the state space, we find that the resulting controls are ro- 
bust with regard to several non-ideal features such as thermal 
fluctuations, decoherence, leakage of the local magnetic field 
and existence of uncertainty in the couplings of the system. 



The paper is organized as follows: in Section II, we briefly 
introduce the model and its invariant excitation subspace char- 
acterized by a definite value of the total excitations. In Section 
III, we will formulate the entanglement generation task into an 
optimization problem, and numerically find the optimal con- 
trol pulses. In Section IV, from controllability point of view 
we provides some analysis trying to understand why such pro- 
posed optimization process can generate the required control. 
In Section V, we add some more practical conditions to the 
problem and discuss how these extra conditions will affect the 
control results. The results are summarized in the section VI. 



II. SYSTEM AND CONTROL MODEL 

We consider an isotropic Heisenberg chain of N spin-^ par- 
ticles with the nearest-neighbor interaction: 

), (1) 

n=l 

where the Pauli operators are defined as usual 



J is the coupling between adjacent spins, and is positive for 
an antiferromagnetic chain. Let |0) and |1) be the basis vec- 
tors of a single spin, representing the spin-down and spin-up 
states. The 2 N -dimensional Hilbert space T~L of the spin chain 
is spanned by tensor products of the TV single spin basis vec- 
tors, the computational basis states, for which we shall use 
the shorthand notation |0 . . . 0) = jO}®^ etc. For each basis 
state the number of excitations is simply the number of | l)s in 
that basis. For instance |11 . . . 1) has N excitations. Setting 

S z = J2n=i Zn> it is eas y to verify that [H , S z ) = 0, hence 
Hq is invariant on every excitation subspaces. Accordingly, H 
can be written as a direct sum of N + 1 excitation subspaces 
Tin (n — 0, . . . , N) whose basis vectors are those having n 
excitations. The dimension of H n is (^). For chains with an 
even number of spins the system ([T| has a unique ground state 
in the y -excitation subspace "Hjv/2- F° r odd-length chains the 
ground state is doubly degenerate with one ground state in the 
H(n-i)/2 and W.(n+i)/2 subspace each. 

Next, we impose a local magnetic field along the z- 
direction on the left end spin, with the resulting Hamiltonian: 

H c = B(t)Hi = B(t)Z 1 (2) 

where B(t) is the magnitude of H c , which can be varied 
over time. For spin chains comprised of coupled quantum 
dots, such an interaction could be achieved using local volt- 
age gates. This is the Hamiltonian we shall control for en- 
tanglement generation. The direction of field is not essen- 
tial due to the symmetry of the Heisenberg Hamiltonian. The 
main idea is to demonstrate that perfect end-to-end entangle- 
ment can be generated from only a local magnetic field along 
a given direction. This requirement is experimentally simpler 
than other proposals where independent control fields along 
different directions are required. Of course, the more control- 
lable degrees of freedom, the more types control tasks can be 
completed, but a single unidirectional local field is sufficient 
for entanglement generation, and indeed controllability on the 
iV/2-excitation subspace as we shall see. 

The total Hamiltonian of the spin chain thus becomes 

H(t) = H S + H C = J Ho + B(t)H v (3) 

Since Hi is invariant on every excitation subspace Hm, H(t) 
is also invariant on these subspaces. As a result the number 
of excitations of a given initial state \ip(t)) does not change 
under the evolution 
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i.e., if |^(0)) is in the AT-excitation subspace then \ip(t)} will 
be in Hm for all t. Throughout this paper we choose units 
such that H = 1, so J and B(t) has the dimension of fre- 
quency. For NMR experiments J is typically a few hundred 
Hz, and B{t) is around 50 kHz for liquid-state NMR and a 
few hundred kHz for solid-state NMR 1 18 1. 



III. OPTIMAL CONTROL DESIGN 

Our question is whether we can find an appropriate control 
pulse B(t) in the model ^ to generate a maximally entan- 
gled state between the two end spins at a certain final time 
t = tf. We shall first assume that the system is initialized in 
the ground state \Go) of Ho and thus in a pure state. In prin- 
ciple, this can be done by simply cooling the systems and this 
assumption will later be relaxed. As mentioned above, for 
even chains the ground state of Ho is unique and lies in the 
subspace ofH.N/2 so that cooling the system is enough to pre- 
pare it. For chains of odd length the degeneracy can be lifted 
by applying a small global magnetic field during the cooling 
process and the system can be prepared into the ground state 
in the subspace T-L(n+i)/2- In the following, for simplicity, we 
will call Hn/2 for even chain and H(n+i)/2 f° r odd chain as 
the largest excitation subspace, and we assume the system is 
initially prepared as the ground state \Gq) in the largest exci- 
tation subspace. 

The control objective is to find a time-varying control field 
B(t), t £ [0, tj], such that the final state at time tf takes the 
form 

|V(t/)> = |^-)l,JV® |^>re, (5) 

where |t/r) = (|01) - |10))/V2 is the singlet state between 
the first and the last spin, and \rp) Te is the state of the remain- 
der, i.e., spins n = 2, 3, N — 1, on which we do not impose 
any constraint. 

One of the methods to generate the control pulse B(t) is 
through optimization, where the target state \ip(tf)) is the 
maximal (or minimal) point of a certain objective function 
K[B(t)]. A good candidate in our case is 

K [B{t)] = Mt/MMi/)), (6) 
where the operator A is the projection onto the target subspace 

A=\i/j-)i !N (ip-\®I r<i) (7) 

and \psi~)i : N{4>~\ is the density operator for the singlet state 
between spins 1 and N, which we desired to prepare, and 7 re 
is the identity operator acting on the rest part of the chain. 
Eq. (|6]l shows that this function is the fidelity of the density 
matrix of the end spins (t) with respect to the singlet state 
at time tf. Consequently, this function is always positive and 
assumes its maximum for states of the form where the two 
end spins form a singlet. Hence, our main goal is to drive this 
system to its maximum by varying B(t) over the time interval 
[0, t f]. For simplicity, we will vary B(t) in the function space 
of the piecewise-constant functions. 



The optimization procedure is as follows. We divide the 
time [0,tf] into p equally spaced sub-intervals, At = tf/p. 
On each time interval the control field B(t) takes the constant 
value B m (m = 1, 2, ..,p) satisfying \B m \ < Bm, where Bm 
is the maximum amplitude of the control field that we can gen- 
erate experimentally. Assuming unitary evolution governed 
by the Schrodinger equation ([3]) the final state subject to this 
piece-wise constant control is simply 

\iP(t f )) = U(B p )---U(B 1 )\G }, (8) 

where 

U(B m ) = e-WWOW ( 9 ) 

and the fidelity is simply a function of p variables, K[B(t)] — 
K(Bi , . . . , Bp). This is a standard optimization problem in a 
finite-dimensional parameter space, which can be solving us- 
ing various methods. We use a quasi-Newton method devel- 
oped by Boyden, Fletcher, Goldfarb and Shanno |20) (BFGS), 
which is an improvement over the standard GRAPE algo- 
rithm 1 19 1 widely used for quantum control problems. 

The algorithm is summarized as follows: (1) Choose an 
initial trial function B(t) = flW(t), which could be constant 
B^- 1 ' (t) = 1 or some other random choice. (2) Iteratively gen- 
erate B^ m+ ^ (t) from B^ (t) by applying the quasi-Newton 
method, which guarantees K[B m+1 (t)] > K[B m {t)]. The 
algorithm is terminated if the change of the fidelity is smaller 
than a given threshold, or we cannot find a search direction 
to increase the fidelity further. The method guarantees that 
the function sequence B^ m \t) converges to a local maximum 
B m ax{t) of the function K[B(t)}, and constraints, e.g., on the 
field amplitudes, can be explicitly taken into account. In ap- 
plications such as NMR where the field amplitudes B{t) can 
be large compared to J (at least 10 2 — 10 3 J) it is usually not 
necessary to impose explicit amplitude constraints — at least 
our simulation results for this problem suggest that the am- 
plitudes of the optimized fields are usually well below 100 J 
if we start with an initial field B^-'{t) close to the zero field 
B = 0, as is to be expected for a local optimization method. 
Therefore in most of the following calculations no amplitude 
constraints have been imposed, except where otherwise stated. 
In general amplitude constraints increase the amount of time 
required to reach the objective but the imposition of such con- 
straints can be beneficial by increasing the robustness of the 
optimal pulses obtained with regard to noise, field leakage, 
etc. Beside the fidelity K[B(t)] of the final state with respect 
to the desired singlet state, we also use the entanglement mea- 
sured by the concurrence [21 1 to quantify the quality of our 
procedure, where the entanglement C is maximal (C = 1) if 
we successfully generate a perfect singlet. 

Using this optimal control framework we computed opti- 
mal pulses for spin chains of length TV = 4 to 20. An example 
of an optimal control field and the corresponding evolution 
of the system for N = 4 is shown in Fig. [T] The control 
steers the ground state of Ho, which is a superposition of all 
six computational basis states in the N/2 = 2 excitation sub- 
space Ti.2, to a +l-eigenstate of the observable A in the same 
subspace. For N — 4 the +l-eigenspace of A restricted to H2 
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FIG. 1 : (Color online) (a) Example of optimal field B z (t) for N = 4 
and (b) corresponding evolution of the populations of the compu- 
tational basis states in the 2-excitation subspace and the population 
p(t) of the state \ip) = |[|0011) - 0101) - 1 1010) + 1100>] (thick 
solid line) and (c) concurrence Ci4(t) between spins 1 and 4. The 
evolution shows that the pulse transfers virtually all population to the 
state j ip) , which exhibits perfect entanglement between the end spins, 
as demonstrated by the concurrence of the reduced density matrix. 
(Pulse shown obtained with amplitude constraint \B(t)\/J < 1.) 



is spanned by {^(|0011) - |1010», ^(|0101) - |H00))} 
and the control must steer the initial state to some state in 
this subspace. In the example shown in Fig. [T] for instance, 
virtually all the population is transferred to the eigenstate 
\tp) = |(|0011) - |1010) - |0101) + |1100))} at the final 
time, although the population of the target state and the con- 
currence between the end spins do not increase monotonically 
unlike for some local optimization approaches El . The par- 
tial trace of the final state over the inner spins (up to three 
digits accuracy) is 



Pi4 = Tr 2i3 (|^)(V'|) 
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which is an almost perfect Bell state with concurrence 
0.99998, and we can achieve fidelities and concurrences al- 
most arbitrarily close to 1 by relaxing the constraints on the 
field amplitudes, pulse duration or time resolution of the fields 
slightly. In this sense we claim that this method can generate 
perfect entanglement. 

Although the dimension of the ^-excitation subspace in- 
creases from (2) = 6 for N = 4 to (*°) = 184, 756 for 
N = 20, rendering the optimization problem considerably 
more challenging and time-consuming, we are still able to 
find optimal optimal pulses for which K[B] assumes values 
very close to 1 in relatively short time, suggesting that our 
optimization method for generating end-to-end entanglement 
is still effective for relatively long spin chains. For a given 
chain there is always a lower bound Tq for the evolution time 



FIG. 2: (Color online) Lower bound on the time Tc required for gen- 
erating entanglement in terms of N 2 , based on numerical optimiza- 
tion results for chains of lengths N = 4 to 15 without constraints. 



tf, however, below which we can never achieve unit fidelity 
no matter how many time steps or iterations are used. Unfor- 
tunately, analytical expressions for these bounds are known 
only for very simple systems. We therefore numerically ex- 
plore the search space to find approximate lower bounds Tc 
as a function of the chain length N. The results, shown in 
Fig.|2]suggest a quadratic increase with respect to chain length 
N, similar to what was reported in [10] for implementing the 
universal quantum computation by controlling two qubits of 
an XY chain, rather than a linear increase as in some other 
entanglement generating schemes [3, 4| where entanglement 
propagates from one end through the chain. The nonlinear 
dependence on the length of the chain is not surprising as 
the mechanism for increasing the entanglement between the 
end spins in our scheme involves dynamic population transfer 
from the ground state of Ho to an eigenstate of the observable 
A. We note that the control itself cannot create entanglement 
directly as the control Hamiltonian Hi is completely local. 
It can only alter the entanglement between the end spins by 
taking advantage of the existing interactions between all spins 
and the resulting non-trivial dynamics. 



IV. CONTROLLABILITY & REACHABILITY 
CONSIDERATIONS 

Considering the restrictions imposed on the control, it is 
rather surprising that we are able to find optimal control solu- 
tions at all, and it is worthwhile to investigate whether there al- 
ways exist controls that achieve perfect entanglement, at least 
in principle. In order to answer these questions, we need to 
review the concepts of reachability and controllability. Ulti- 
mately, the question we would like to answer is what is the set 
of states that is reachable from the ground state for a system 
governed by the evolution Q. We already know that there 
are some constraints. For a Hamiltonian system the spectrum 
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of the density operator is invariant. For pure initial states this 
simply means that the set of reachable states is a subset of 
the pure states. We also know that there are additional limita- 
tions as both Hq and Hi in our case commute with total ex- 
citation operator S z , and hence the Hilbert space decomposes 
into excitation subspaces on which the dynamics is invariant. 
Thus if the initial state is a pure state on the ^ excitation sub- 
space then the set of reachable states is a subset of the pure 
states on this excitation subspace. However, the key question 
is whether the target state Q in this excitation subspace is 
reachable from the initial state \Go). Unfortunately, this is a 
very difficult question to answer directly in general, except for 
controllable systems. We must be careful here as there are dif- 
ferent degrees of controllability [23| but if we can show, for 
instance, that the unitary evolution of the system is control- 
lable, then pure-state controllability follows. 

The Schrodinger equation Q can be rewritten into an equa- 
tion in terms of the unitary process U(t) as the variable 

U{t) = -t(JH + B{t)H x )U{t) (10) 

where U(Q) = I and the solution for a given control B(t) is 
denoted as U(t, B(t)). In this context, the reachable set TZ is 
defined as the set of unitary matrices that can be achieved by 
the dynamics U (t, B) for some (admissible) control B(t), i.e., 
U € TZ if and only if there exists an admissible control B(t) 
such that U(tf, B(t)) = U and we have the following [24 1 : 

Theorem 1. Let C be the Lie algebra generated by 
span{— iHq, —iHi}, often called the dynamical Lie algebra. 
Then TZ = e c . 

Thus, the reachable set depends on the dynamical Lie al- 
gebra generated by the Hamiltonian. If C = u(2 JV ) or 
C = su(2 JV ), then TZ = V(2 N ) or TZ = SV(2 N ) and the 
system is controllable in the sense that we can implement any 
unitary process up to a global phase and steer any density op- 
erator to any other density operator with the same spectrum. 
Back to our spin chain problem, for H = JH + B(t)H\ 
with Hq and Hi as in Eqs ([T]l and if the system is control- 
lable on H.n/2 for N even or H(n±i)/2 f° r N odd, then the 
existence of the control to generate the required singlet state 
is guaranteed. A straightforward way to check controllability 
is by calculating the dimension of the dynamical Lie algebra. 
Calculations for small chains of different length suggest full 
controllability on the relevant excitation subspaces for chains 
of both even and odd length. However, such calculations are 
laborious and impractical for long chains. Similarly, prov- 
ing that the dimension of the dynamical Lie algebra is M 2 , 
where M is the dimension of the TIn/2 ° r ^(jv-i)/2 excita- 
tion subspace, for any N is a difficult mathematical problem. 
However, for even-length isotropic Heisenberg chains, con- 
trollability can be easily verified using the following useful 
result Il25l : 

Theorem 2. If Hq is effectively strongly regular and H\ is 
connected, then the system is controllable. 

Here Hq is strongly regular if all the transition frequen- 
cies uike = efe — ee, where are the eigenvalues of Hq, 



are distinct. Choosing a basis such that Hq is diagonal, 
Hq = diag(ei, e2, . . . , e„), let Hi = (fry) be the matrix 
representation with respect to this basis. Hi is connected if 
the graph represented by Hi is a connected graph. To be 
more specific, we say the indices j and k are connected if 
bjk 7^ 0, and if every two indices are connected through sev- 
eral nonzero 6,-fc's, then Hi is called connected. Moreover, 
Hq is called effectively strongly regular, if uijk, j 7^ k, are 
nonzero and distinct whenever bjk 7^ 0. Using these defi- 
nitions and the previous theorem we can verify on a case by 
case basis that even-length Heisenberg chains are controllable 
on H N/2 - 

Example 1. For N = 4 we have dim(%) = 6, Hq and Hi 
take the following forms in the eigenbasis of Hq: 



Hq = 



Hi = 
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which shows that Hq is effectively strongly regular and Hi 
is connected. Therefore, the system is controllable on the 
H2 excitation subspace, and by virtue of controllability any 
pure state in tlie Tii excitation subspace is reachable from the 
ground state \Gq), including our target Bell state ([5]). 

We can similarly verify that the systems governed by ([3]) 
satisfy the criteria for controllability above on Hn/2 f° r N 
even ranging from 4 to 16. For longer chains, such verifica- 
tion will become inefficient again due to the evaluations of 
eigenvalues and eigenvectors for large matrices. For odd N, 
Hq ceases to be effectively strongly regular on the (N ± l)/2 
excitation subspace for N > 3, and the above theorem does 
not apply. However, this does not imply non-controllability 
as the above controllability criteria is only sufficient but not 
necessary. In fact direct calculation of the dimension of the 
dynamical Lie algebra for odd chains suggests that the sys- 
tem is still controllable. For instance, for N = 5 we obtain 
dim 7^2 = (2) = 10 and direct computation of the Lie alge- 
bra shows that it has dimension 100 = 10 2 , suggesting that 
the Lie algebra in this case is u(10), which implies full con- 
trollability on the subspace. 

Although controllability implies reachability and the ex- 
istence of a control that steers the initial state to the target 
state, it should be noted that it does not answer the question 
of the minimal time required to steer the system to the tar- 
get state. Furthermore, in practice the finite time resolution 
of the field and bounds on the field strengths impose addi- 
tional constraints, i.e., even if there exists an unbounded con- 
trol that steers the initial state to the target state in time tf, 
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there may not be a piecewise constant solution of the form 
(B m ) = (Bi, . . . , Bp) if p is fixed and \B m \ < Bm for some 
fixed bound Bm- Finally, most optimization techniques in- 
cluding the quasi-Newton-method used in the previous section 
are only guaranteed to find some extremum of the objective 
function, which need not necessarily correspond to a global 
optimum. 

For more general spin chain models, for example, the XXZ 
chain, we can still use the above optimization method to gen- 
erate perfect entanglement. For different length of chains we 
have calculated, the system is always found to be controllable 
on the subspace in question. More interestingly, for XX chain, 
such optimization method can also generate perfect entangle- 
ment from the ground state, but the system is found to be non- 
controllable on that subspace. This gives an example where 
the perfect end-to-end entangled state is reachable from the 
ground state, but controllability does not hold, a strong evi- 
dence that reachability is weaker than controllability. 



V. CONTROL PERFORMANCE UNDER REALISTIC 
CONDITIONS 

In the previous sections we have illustrated that a single 
magnetic field in the z-direction acting on a single spin suf- 
fices to generate a perfect entanglement between two ends of 
a Heisenberg chain. However, the optimal control design is 
based on four ideal assumptions, namely that 

1. The system can be prepared in its ground state; 

2. The control field can locally address a single spin with- 
out affecting the others; 

3. The system is homogeneous, i.e., all couplings are equal 
and the coupling strength is known; 

4. The system is isolated from its environment. 

In this section we consider how the effectiveness of the opti- 
mal controls is affected when these assumptions are relaxed. 



A. Thermal Initial State 

In an ideal setup we would be able to cool the system to its 
ground state |Go) but in practice this is difficult, and thus the 
initial state is likely to be a thermal ensemble 



P(0) 



-H a /k B T 



Tr[e- ff o/fe B T] ' 



(11) 



where T is temperature and ks is the Boltzmann constant, 
which contains some contribution of the excited states, in- 
cluding states outside the subspace Hn/2- To assess the effec- 
tiveness of the optimal control pulse generated for the ground 
state in this case, we evolve the initial state ( fTTj ) under the ac- 
tion of the total Hamiltonian H (t) given in Eq. ^ subject to 




k B T/J 



FIG. 3: (Color online) Final concurrence between end spins vs 
kBT/J for chains of different lengths if the pulse optimized for en- 
tanglement generation from the ground state is used. For large fcsT 
the concurrence approaches zero and the pulse is ineffective, but for 
fcsT/J <C 1 (low temperature limit) the concurrence remains close 
to unity. 



the pulse B max (t) optimized for the ground state, 



p(t f ) = exp H 



H(t)dt 



p(0) exp_, 



! t 



H(t) dt 



(12) 

where exp + indicates positive time ordering, and compute the 
density matrix for the end spins pi n (tf ) by tracing out the in- 
terior spins. Fig. [3] shows the resulting entanglement as mea- 
sured by the concurrence as a function of fc^T. For suffi- 
ciently low temperatures T the optimal control field still gen- 
erates good entanglement. For example, for a chain of length 
N = 10 entanglement survives up to fcsT/J ~ 1.2. More- 
over, Fig. [3] shows that for low temperatures the concurrence 
is rather flat, implying that up to a certain temperature the con- 
trol performance is minimally affected by thermalization. The 
width of the plateau is related to the energy gap between the 
ground state and the first excited state. For thermal energies 
kgT less than this energy gap the population of the excited 
states remains very small, and only when fc^T exceeds this 
limit other eigenstates become populated. Since the energy 
gap between the ground state and first excited states decreases 
with chain length, the width of the plateau also decreases. 

These results are consistent with expectations. For low tem- 
peratures the ground state of the system dominates the mix- 
ture of the thermal state (Hi, and the optimal control pulse 
transfers this to a +l-eigenstate of A that has unit concur- 
rence. Most of the other states of the initial ensemble will be 
mapped to 0-eigenstates of A and not contribute to the entan- 
glement between the end spins, but if the initial populations 
of the excited states are small, the pulse will remain mostly 
effective. The results can be improved greatly if we calculate 
the optimal pulse for the initial thermal state, rather than us- 



7 




E 0.7 



FIG. 4: (Color online) (a) Optimal control pulse B z (t) obtained 
by optimizing Tr[Ap(tf)] for an initial thermal ensemble with 
k B T/J = 2 for N = 4. (b) If if this pulse is applied to ini- 
tial ensembles p(0) with fcsT/J ranging from to 2, the fidelity 
Tr[Ap(t f)] achieved (squares) assumes its theoretical upper bound 
(solid blue line), showing that this field is optimal in terms of maxi- 
mizing Tr[Ap(tf)] for initial ensembles over a wide range of temper- 
atures. Accordingly, the corresponding entanglement as measured 
by the concurrence of the reduced density matrix p\$(tf) for the 
ensemble-optimized pulse (dashed red line with diamonds) is sub- 
stantially higher than that for a ground state optimized pulse (green 
dash-dot line with circles). 



where we have ordered the spectrum and the correspond- 
ing eigenstates in an non-increasing way for both p(tf) and 
A. Different thermal ensembles at different temperature 
have different spectra w n , m - However, to realize the upper 
bound ( 13 I, any optimal pulse for a given thermal ensem- 
ble must map all the eigenstates |e„. m ) to the eigenstates 
of A in the same order, independent of the actual value of 
w ntm . Therefore, the optimal pulse for a thermal ensemble 
with 7i > is also an optimal pulse for any thermal ensem- 
ble with T2 < T\, and in particular, it is an optimal pulse for 
the initial ground state. But the inverse is not true: the opti- 
mal pulse obtained if we start in the ground state is in general 
not optimal for an initial thermal ensemble. This is because 
the optimal field starting in the ground state is only guaran- 
teed to map the final state onto one of the +l-eigenstates of 
A, but it generally will not map the correct excited states onto 
other +l-eigenstates of A, and therefore it is not an optimal 
pulse for an initial thermal ensemble. This is illustrated in 
Fig. |4j which shows that the fidelity that can be achieved if a 
ground-state optimized pulse is applied to a thermal ensemble 
with T > is strictly smaller than the optimal fidelity that is 
achievable with a pulse optimized for this thermal ensemble. 
The concurrence exhibits a similar behavior. For N = 4 the 
final concurrence between the end spins for ksT/J = 1 is 
0.9762 vs less than 0.7 for the ground-state optimized pulse 
in Fig. [3] 



B. Local control with Gaussian leakage 



ing the optimal control for the ground state, i.e., maximizing 
the fidelity Tr[Ap(fy)] starting with a thermal ensemble p(0), 
although there are upper bounds on the maximum achievable 
fidelity in this case. Let p(0) = ®^ p {n} (0) where p^> (0) is 
the initial state restricted to the nth excitation subspace, w n , m 
be the populations of the initial ensemble on the nth excitation 
subspace with w n , m > w rl/m +i > for rn = 1, . . . , dimH n , 
and let d n be the number of +1 -eigenvalues of A on the nth 
excitation subspace. Then the constraint of unitary evolution 
combined with the fact that the dynamics is invariant on each 
excitation subspace imposes the following upper bound on the 
maximum achievable fidelity |26l 

N d n 

Tr[Ap{t f )] < £ Wn,m- (13) 

n— m— 1 

Fig. |4] shows that we can effectively realize this upper bound 
for a wide range of fcgT. In fact, the figure shows that it 
is not necessary to know the initial temperature. The pulse 
that achieves the maximum fidelity possible for ksT/ J = 2 
is an optimal solution for all temperatures in that it achieves 
the maximum possible value of Tr[^4p(fy)] for every initial 
thermal ensemble p(0). This can be explained as follows. De- 
noting the eigenstates of Hq by \e njn ) for n = 0, . . . , N and 
m = 1, . . . , dirnHn, the initial state can be written as 

JV dim«„ 

^ ^ ^ ^ ^n.m |^n,m) (^n,m | ; (14) 
n— m— 1 



Another central assumption of our control scheme is that 
the control field is confined and interacts only with the first 
spin. In practice such perfect confinement is difficult to 
achieve and any local magnetic field will also affects neigh- 
boring spins, although the effect will usually decay with in- 
creasing distance from the control spin, the first spin in our 
case, either exponentially or following a power law. 

To assess the effect of leakage on the performance of the 
control scheme we take the pulses optimized under ideal as- 
sumptions and consider what happens if the field decays expo- 
nentially with distance thus affecting not only the target spin 
but also neighboring spins, effectively changing the control 
Hamiltonian to 

N 

H c = B(t)H 1 = B(t)Y,e-< k - m Z k (15) 

k=l 

where £ is a parameter quantifying the field leakage with 
£ = corresponding to the ideal case where the field is com- 
pletely localized and affects only the first spin. Leakage in- 
creases with £ and for £ —> oo the control field is uniform, 
resulting merely in a dynamic global energy shift or global 
phase. Fig.|5]shows the effect of leakage on the entanglement 
between the end spins at the final time as a function of the 
parameter £. As expected, entanglement decays with increas- 
ing £ but even for 10% leakage we still have a good entan- 
glement. Furthermore, if the leakage parameter £ is approxi- 
mately known, we can again incorporate this information into 
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FIG. 5: (Color online) Entanglement as a function of leakage £ for 
different chain length (N = 6 black triangles, N — 8 blue squares, 
N = 10 red circles). Solid lines indicate the concurrence depen- 
dence for pulses optimized assuming ideal Hi and dashed lines in- 
dicate concurrence behavior for pulses optimized with leakage taken 
into account, i.e., Hi. For £ small the optimal pulses obtained as- 
suming no leakage still achieve high concurrence but if £ can be esti- 
mated, the entanglement generated can be substantially improved by 
taking leakage into account in the optimization. The inset shows that 
the optimal pulses in the presence of leakage can differ substantially 
from the optimal pulses in the absence of leakage. Note that opti- 
mal control solutions are not unique and the robustness of individual 
solutions with regard to leakage varies. 



the optimization, replacing Hi by Hi and then doing the op- 
timization. This substantially improves the results, allowing 
us to achieve consistently high entanglement for a wide range 
of £, although for very large £ Hi almost commutes with Ho, 
rendering the control increasingly ineffective (see Fig. [3J. 



C. Random Couplings 

Another assumption in our model-based optimal control 
scheme is that the system Hamiltonian Ho is homogeneous 
and known. In practice, there is no guarantee that the Hamil- 
tonian of the system will be this perfect. In particular the 
J-couplings between neighboring spins are likely to be sub- 
ject to random variations and the precise actual values of the 
couplings may not be known, leading to an uncertainly in the 
system Hamiltonian. It is therefore very important to con- 
sider how small random variations of the coupling strengths 
are likely to affect the effectiveness of the control scheme. To 
this end, let us assume that the coupling between neighboring 
sites are random and vary around an average value J, leading 
to an actual system Hamiltonian of the form 

JV-l 

H s = JHq = J (l+e n )(X n X n+ i+Y n Y n+ i+Z n Z n+ i) , 

71=1 

(16) 




0.03 0.06 0.09 0.12 0.15 
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FIG. 6: (Color online) Average entanglement, (Cijv), in terms of 
disorder parameter a for different lengths. Each point is the average 
entanglement over 100 different randomly perturbed Hamiltonians. 
A small amount of uncertainty in the coupling constants hardly af- 
fects the entanglement achieved but for larger variations the entan- 
glement drops significantly. The average entanglement achieved for 
disordered systems when the actual couplings are used in the opti- 
mization (solid black line with triangles for N — 6) is very close to 
one even for high disordered systems, showing that inhomogeneity is 
not a limiting factor provided the actual couplings can be estimated 
with sufficient accuracy. 



where e„ E [—a, +a] is a uniformly distributed random vari- 
able with mean 0. Assuming that the precise actual couplings 
are not known, the best we can do is to study the average ef- 
fect of random variations of the coupling strengths of a cer- 
tain magnitude a on the entanglement generated. To this 
end, we evolve the ground state of the (randomly) perturbed 
Hamiltonian ( [To*} according to the Schrodinger equation with 
H = JHq + B(t)Hi with B(t) being the pulse optimized 
for the ideal Hamiltonian Hq and its ground state \Gq). The 
resulting entanglement between the end spins will usually be 
less than in the ideal case. For comparison we generate 100 
different random Hamiltonians Ho for each value of a, com- 
pute the entanglement attained between the end spins as quan- 
tified by the concurrence Cijv if we apply the optimal pulse 
for the ideal case, and take the average (Cin)- 

Fig. [6] shows the average entanglement as a function of a. 
Although the average entanglement decreases with a, it re- 
mains high for variations up to a = 0.01J, showing that 
a small amount of uncertainty can be tolerated. Larger un- 
certainties in the coupling strength can be addressed by ex- 
perimental characterization using, e.g., spectral analysis or 
Bayesian estimation techniques [27], and the figure shows that 
it is highly beneficial to experimentally characterize the cou- 
plings to reasonable accuracy, or alternatively, to use adaptive 
closed-loop experiments (111 to refine the results as disorder 
or inhomogeneity is not a limiting factor in the optimization 
per se. If the coupling constants used in the optimization are 
sufficiently close to the actual couplings, near perfect entan- 
glement can be achieved even for inhomogeneous systems. 
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Another interesting feature of Fig. [6] is the fact that longer 
chains are less sensitive to the random couplings, and thus do 
not need to be characterized with the same accuracy necessary 
for shorter chains to obtain the same amount of entanglement. 
This can be explained in terms of fundamental mathematical 
properties and the central limit theorem [28 1, which asserts 
that when iV identical random variables are added, the dis- 
tribution of the result is Gaussian and its mean and standard 
deviation are normalized by N and 1 / yN respectively. This 
means that the randomness of the system Hamiltonian ( [To} is 
effectively suppressed by a factor of \/N with increasing N. 

D. Decoherence 

Another source of errors is from decoherence and popula- 
tion relaxation induced by incoherent interactions of the sys- 
tem with the environment. On the short time scales considered 
here we expect decoherence due to pure dephasing to be the 
dominant effect and the system evolution to be governed by a 
master equation of Lindblad form 

p=-i[H(t),p]+£(p). (17) 

To assess how such environmentally induced phase relaxation 
affects the performance of the optimal control scheme, we 
compute the evolution of the system according to (JTTJi starting 
with the initial state po = \Go)(Gq\, subject to the Hamil- 
tonian H(t) = JH + B(t)Hx as in (TJl, where B(t) is the 
optimal pulse for the ideal system, and again calculate the 
concurrence between the end spins of the final (mixed) state 
P(tf)- 

Assuming that the main source of phase relaxation is from 
the incoherent interactions with the control and measurement 
apparatus, we consider two types of phase relaxation de- 
scribed by the Lindblad operators 

A(p) = -7{2p - ZipZi - Z NP Z N }, (18a) 

N 

-CjW = -7E^-W (18b) 

In the first scheme dephasing affects only the end spins, which 
are interacting with the probes, while in the second scheme we 
assume that all spins are affected by dephasing due to imper- 
fect shielding. The resulting entanglement between the end 
spins in terms of noise strength 7 in both dephasing cases is 
shown in Fig. [Tjfor a chains of length N = 6 and A" = 10. 
Decoherence reduces the entanglement but when the noise 
strength is small enough we are still able to generate entan- 
glement with the pulse optimized for the ground state. When 
more spins are affected by decoherence as in the second sce- 
nario, the entanglement decays faster, and the effect of de- 
coherence increases with chain length, as expected. For the 
scheme to be effective if all spins are subject to dephasing 
for a chain of length N = 10, the relative decoherence rate 
7/ J should be less than 0.1%, which is consistent with usual 
quantum information processing requirements. The results 
can generally be improved somewhat if decoherence is taken 
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FIG. 7: (Color online) Entanglement as a function of dephasing rates 
7/ J for chains of length N = 6 and N = 10 for dephasing on 
the end spins only (a) and dephasing of all qubits (b) as defined in 
Eq. j 1 8ft for a pulse optimized for the Hamiltonian system, and for 
TV = 6 concurrences achieved with decoherence is taken into ac- 
count in the optimization (*). 



into account in the optimization process but there are limita- 
tions as coherent control cannot compensate for irreversible 
effects caused by the Lindblad operators. For 7/ J <C 1 we 
find that the pulses obtained by optimization over the Lindblad 
dynamics generate almost the same amount of entanglement 
as those computed assuming Hamiltonian dynamics, and for 
longer chains improvements seem to be even less although it 
must noted that there are many solutions to the optimal control 
problem for Hamiltonian systems that achieve near unit con- 
currence, and the robustness with regard to decoherence of 
these pulses may vary. Thus a more extensive global search 
of the parameter space may yield better results. 



VI. CONCLUSION 

We have studied control of anti-ferromagnetic isotropic 
Heisenberg chains using a single control field such as a mag- 
netic field along a given direction acting on a single site. The 
intrinsic isotropic Heisenberg Hamiltonian commutes with the 
total spin operator, implying that the dynamics decomposes 
into N + 1 independent excitation subspaces. Although the 
control is not sufficient for full controllability of the system 
as the resulting control Hamiltonian commutes with the total 
spin operator, and hence each excitation subspace is invariant 
under the dynamics regardless of the control field applied, we 
verified that the system is controllable in the largest excitation 
subspace, which guarantees the existence of a control driv- 
ing the ground state to a perfect end-to-end entangled state in 
the subspace. Such control can be generated by transforming 
the control problem to an optimization problem and deriving 
the optimal control through standard optimization procedure. 
More importantly, it was shown that the resulting optimal con- 
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trol solutions are quite robust with regard to multiple sources 
of imperfections, including imperfect initial state preparation, 
incomplete confinement of the local control fields, uncertainly 
about the couplings between adjacent spins and decoherence 
due to environmental influences. 

Furthermore, the quality of the entanglement achieved can 
be substantially improved by taking the imperfections into 
account in the optimization. If we can estimate the amount 
of leakage of the control field or the actual values of the 
couplings between adjacent spins using experimental system 
identification techniques, then such information can be incor- 
porated in the optimization to generate greatly improved op- 
timal pulses. Furthermore, we found that a control pulse that 
maximizes the entanglement between the end spins starting 
with a high-temperature thermal ensemble is automatically 
optimal for lower-temperature initial ensembles, which means 
that the optimal control scheme is effectively independent of 
the initial state, provided it is a thermal ensemble. Another in- 
teresting observation is that simulations strongly suggest that 



for our problem controllability on the largest excitation sub- 
space holds not only for the isotropic Heisenberg chain but 
also for a general XXZ model, which deserves future investi- 
gation. 
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